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Abstract. We consider the asymptotic profiles of the nonlinear parabolic flows 

(e") t = Au + Ae" 

to show the geometric properties of the following elliptic nonlinear eigenvalue problems known as a 
Gelfand's problem: 

A<p + Ae'' = 0, cp > in Q 
<p = on Q 

posed in a strictly convex domain OcR". In this work, we show that there is a strictly increasing 
function f(s) such that / _1 (<p(x)) is convex for < A < A", i.e., we prove that level set of cp is convex. 
Moreover, we also present the boundary condition of cp which guarantee the /-convexity of solution 
cp. 

1. Introduction 

We will investigate the geometric properties of parabolic flows and derive related geometric 
properties for the asymptotic limits of such evolutions. More precisely, we consider the nonnega- 
tive solutions u(x, t) of the following equation 

(1.1) (e'% = au + Ae 11 

posed on a strictly convex and bounded domain Q with zero boundary condition 

(1.2) u = on dQ., 
and initial data 

(1.3) u(x, 0) = u Q {x) > 

In the limit, these flows converge to solutions of the well-known gelfand's problems 

fA<p(x) = -Ae'PW inQ 
{(p(x) = on dCl. 

By this relation, it is natural to expect that the solution of the parabolic flow above have a lot in 
common with those of the Gelfand's problems, ( |GPa| >. The aim of the paper is to provide geometric 
properties of solutions of the Gelfand's problems by using the prabolic method which is introduced 
by Lee and Vazquez, ILVI . 

Recently, there are a lot of studies of the problem \GPa I because of its wide applications. It arises 



(GPa) 



in many physical models: it describes problems of thermal self -ignition IGel , a ball of isothermal gas 
in gravitational equilibrium proposed by lord Kelvin IChl , the problem of temperature distribution 
in an object heated by the application of a uniform electric current [KG I and Osanger 's vortex model 
for turbulent Euler flows IICLMPI . 
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Let us summerize some known results to {GP^ I. It is well-known that there exists a finite 



positive number A*, called the extremal value, such that Gelfand's problem ( |GPa| ) has at least a 
classical solution which is minimal among all possible positive solutions if < A < A*, while no 
solution exists, even in the weak sence, for A > A*. Let us call the minimal solution cp . The family 

of such solutions depends smoothly and monotonically on A, and in particular 

cp < cp ifA<A'. 
Ha -LA' 

Their limit as A / A' is the external solution cp\- and it can be either classical or singular (i.e. 
unbounded). It is known that cp* e L TO (Q) for every Q if N < 9, while cp* - log (^p ) and A* = 2(n - 2) 
if N > 10 and Q = B\. Brezis and Vazquez |BVj investigated the existence and regularities of 
extremal solutions when they are unbounded. The regulartiy thoery for the Gelfand's problem 
was improved by S. Nedev [Nell who proved that, for general smooth domain Q, cp* is a classical 
solution if n < 3, while cp* e Hq(Q) if n < 5. 

The ultimate goal in this article is to establish the geometric properties of cp for A < A*. Especilly, 
we'd like to show that 

f(cp) - e 5 -a : convex 
From now on, we call it /-convexity shortly. We also refer to the minimal solution cp of the 
Gelfand's problem ( GPa| ) as cp for the rest of the paper. Since the minimal solution cp can be 



obtained as the limit of solution u to (|1.H> -(|1.3 |) as t — > oo, we will concentrate on showing f- 
convexity of u under the assumption that the initial value Uq has the following property 

/(wo) = e"5"° : strictly convex, 

The parabolic approximation method introduced in IILVII relies on the fact that the nonlinear 
elliptic problem ( |GPa| ) can be describe the asymptotic profile of a corresponding parabolic flow 
in a bounded domain and we use that possibility as follows: we select an initial data for the 
parabolic flow having the desired geometric property. Then, the correspoonding solution u of 
will converge eventually to the minimal solution cp as t — > oo. If the evolution preserves the 



/-convexity property under investigation, the result for the problem (|GPj|) will be obtained in the 
limit t — > oo. 

To investigate the /-convexity of solutions to the corresponding parabolic flow, we will split it 
into two steps. The step 1 will be devoted to the study of the /-convexity of solution u on the 
boundary. On dQ, the second derivatives of f(u) can be written in the form 

1 _i„/l 



[/( W )U = 2 S ^ ( 2 M " ~ Uaa ) ' ( - De " M = Ua) ' 



Hence, the geometric properties of solutions on the boundary can be determined by the balances 
between quantities, u a and u aa . However, the difference between them is very subtle in this prob- 
lem. Therefore, there is little room for perturbing quantities. Thus, it is very difficult to investigate 
the geometric properties of solution on the boundary without some boundary condition. In this 
paper, we will focus on the solutions of having the conditions not only {1.2\ , dl -3b but also 

(1.4) G(u, A, Q) = \i\ + A + (n - l)u v H(dQ.) + Ku v > on d£t 

for sufficiently large K > where v is the outer normal vector to dQ, and H(dCl) is the mean 
curvature of Q at the boundary. Here the constant K is related to the shape of boundary <9Q. 

As the second step, we extend the geometric result on the boundary to the interior of Q. Since 
the continuity of the second derivatives of /(«), it is expected that if there are some problems or 
difficulties then they may occur at the region far away from the boundary. Hence, the equation that 
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describe the parabolic flow plays a very important mathematical role in the study of the geometric 
properties of solution in the interior. 

1.1. Outlines. This paper is divided into three parts: In Part 1 (Section 2) we study the convexity 
of solution to the degenerate equation 

(1.5) w t = p(w)a i] \Vw)wij + b(w)\Vwf + c(w)\Vw\ 2 + d{w) 

on a strictly convex domain Q. It is a simple observation that the convexity of solution will 
be strongly effected by the coefficients. Hence, proper conditions need to be imposed to the 
coefficients of ( 11.51 ) for the result. The Part 2 (Section 3) is devoted to the proof for the /-convexity 
of the minimal solution to the Gelfand's problem with boundary condition \1A) . As will be 
mentioned later, the main equation of Gelfand's problem is a special form of dl -5b - Thus, the proof 
will be focused on the /-convexity of solution on the boundary. In Part 3, we will discuss the 
boundary condition \\A) of a solution 

2. Convexity for degenerate equation 

In this section, we will study degenerated equations of the form 

(2.1) w t = p{w)a i \Vw)w ij + b(w)\Vwf + c(w)\Vw\ 2 + d{w), (p > 2) 

on the bounded cylinder Q X [0, oo), where Q is a bounded strictly convex domain in R" with 
smooth boundary. The subindices i, j € {1, • • • , n] denote differentiation with respect to the space 
variables X\, • • • , x n and the summation convention is used. We assume that the coefficient matrix 
(a l i) is strictly positive and all coefficients a l i, b, c and d belong to appropriate C k , (k = 1, ■ ■ ■ ) space 
which will be defined later. The degeneracy of the equation is carried through the function p(w) 
which is assumed to be smooth on Q. 

We assume further that the coefficients of (|2.1|) satisfies the following conditions: 
I.l b(s), c(s) and d(s) are convex, 

1-2 9"-%t = 0. 

2p 

Denoting by L the operator 

Lw = w t - [p{w)a l ^w)Wij + b(w)\Vwf + c(w)\Viv\ 2 + d(w)) , (p > 2) , 
we can now state the main result in this section: 

Lemma 2.1. Let Clbea strictly convex bounded domain in R" with smooth boundary and suppose that the 
coefficients a'i, b, c and d of the operator L are smooth and satisfy the elliptic condition 

a !i U, > c |£| 2 > V£ € R"\{0) 

for some positive constant Co In addition, assume that p is a smooth function on Q and strictly positive in its 
interior. Let wbea smooth solution of (|2.1|) satisfying the conditions conditions 1.1-2. Ifw is convex on the 
parabolic boundary o/Qx [0, oo), i.e., ifmi dinmoo) inf u= i w aa (x, t) > and inf Q rru> a |=i w aa (x,0) > 0, 
then w is convex in the space variable for all t > 0, i.e., infQ X [o /00 ) inf| ea |=i w aa (x, t) > 0. 
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Proof. By direct computation, (|2.1[) implies that the evolution of w a p is given by the equation 
(2.2) 

w A pt =p{w)a ]i {yw)w a f ii ij + p(zv)a^07zv)(w ak Wpij + wp k w ai j) 

+ p'(w)a l i(Vw)(w a Wfsij + WpWaij) + p' (w)w a a l >w k pDjjiv + p'(w)wpd^w ka Wij 

+ p{w)dl{Vw)w kafi Wii + p(w)a'l l (Vw)w ka wipDijiv 

+ p' (w)w a pa'i (Viv)Dijiv + p"(w)w a zvpa' ] (Viv)DjjW + b'\Vwfw a p 

+ [pfc(w)|Vw|P" 2 + 2c{wj\ (Vw a ■ Vwp + Vw ■ Vw a p) 

+ p(p - 2)b(w)w k WiW ka 'Wip |Vwf~ 4 + [pb'(zv )\Vwf~ 2 + 2c' (zf )] (w a Vzv ■ Vzvp + zopVzv ■ Vw a ) 

+ b" (zv)\Vw\ p ~ 2 w a zvp + c' (iv)iv a p\Vw\ 2 + c"(w)w a wp\Vzv\ 2 + d'{w)w a p + d"(w)zv a zvp 

To estimate the minimum of the second derivatives of zv with respect to space variables, we take a 
look at the following quantity, for a positive function ip(t), 



inf inf \waa + eip(t)] = w m {x , f ) + ei/>(f ) 
□,se[0,f] egeRMeeNl L 



xeQ,se[0,f] e^eKM^I 

We need to show that there are a function i/>(f) and £q such that 

w«a(x , fo) + #(fo) > VO < € < £q. 
To get a contradiction, suppose that 

(2.3) w m (x Q) t ) + eil>(t ) = 0. 

Observe that the minimum is taken at (xo, to) with direction a. By the previous lemma, the minimum 
point of iv aa is little way off the boundary. Hence, we can put Xq without loss of generality. On 
the other hand, the parabolic equation w aa contains third order derivatives which is difficult to 
contral with the information of zv aa . Hence we are going to perturb the direction of the derivative 
to create extra terms, keeping the minimum point and minimum zero. 
We now use the function 

Z = w A pr } a ^+e^(t)\r 1 \ 2 

where the modifying functions rf{x) are constructed as follows: at x — we assume that the rf 
satisfy the system 

(2.4) rf. = {c y rf)5pi, rf.. = (cpffybyfyj, 

where the subscripts are space derivatives. Putting also rf{$) - b^p, if follows that 

and 

rf(x) = b^p + C&P + -CaCyXYx? . 

Hence, at x - 0, 

Z ri = w, a pirfif + 2w p c a r] a rf + eip(t)(\t]\ 2 ). 
Z,ij = w A pifrfrf + 2c a w i pijrfrf + 2cpw Ai jrfrf 

+ 2cf a w i p i tfrf + 2c i c a w i pftfrf + 2c a CpW tij rfrf + eip(t) (|r/| 2 ).. . 



(2.5) 
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Zjt - eip'\i]\ 2 =pa ! 'DijZ + pd[ [w ak WpijT] a rf + w^ k w ai jifrf\ 

+ a'i [p'w a - 2pc a ] WpijTfrf + a'i [p'wp - 2pc^] w a ijtfrf 
+ a 1 ' [p"w a wp - 2pc a cp[ Wijtfrf - 2p^cf a w^rfrf - Ipa^qcaW^ifif 

- eipa'i (|j]| 2 ).. + p^winftWijrfrf + p^w^w^w^frf 

+ p'WapOijWijlflf + p'ajwaW^Wijlflf + p'a^WpW^Wijlflf 

+ [pb(iv)\Vwf- 2 + 2c(w)] (Vw ■ VZ) + p\Vwf~ 2 \{b'w a - c a b) Vw ■ Vtoptfrf] 
+ p\Vivf- 2 [(b'wp - cpb) Vzv ■ Vw a ifif] + Uc'w a - c a c) Vw ■ VwpTfrf] 
+ Uc'wp - cpcj Vw ■ Vw a rfrf\ - eip \pbiw)\Vwf~ 2 + 2c(w)] (Via ■ V|r]| 2 ) 
[pb(w)\Vwf~ 2 + 2c{w)\ (Vw a ■ Vivp) if if + p{p - 2)b{w)\Vwf- i w k WiW ka Wiprfrf 
+ [b"\Vwf- 2 + c"\Vw\ 2 + d"\ Wawpfrf + [b'\Vwf + c'\Vw\ 2 + d'] w a prfrf. 
We are now going to choose c a such that 

(2.6) p'w a - 2c a p = 
at (xq, to). Then, by the conditions I and (12.51 1, 

Z, t - eip'\if >pa l iDijZ + pa* \w ak w^rfrf + Wp k w ai jifr]^ 

- IpttkfaWfiTfrf - 2pa % iciC a Wfijtfrf 

- eipa l i (\t]\ 2 ).. + pa'lwij [z k - 2c a w^ k rfrf - eip (l?]l 2 )J + ptfyokaatyWijifif 
+ p'wapdhvifcfrf + p'afaaWtyWijrfrf + p' ^wfa^Wiffrf 

(2.7) + [pb(w)\Vwf- 2 + 2c(w)] (Vw ■ VZ) + p\Vwf~ 2 \{b'w a - c a b) Vw ■ VwpTfrf] 
+ p\Vwf~ 2 [(b'zvp - cpb) Vw ■ Vw a rfrf\ + [{c'w a - c a c) Vw • Vrvptfrf] 

+ [(c'wp - cpc) Vw ■ Vw a ifif\ - eip [pb{w)\Vwf- 2 + 2c{w)\ (Vw • V\tf) 
[pb(w)\Vwf~ 2 + 2c(zv)] faw a ■ Vwp) if if + pip - 2)biw)\Vwf~ i w k wiw ka wipifif 
+ [b'\Vivf + c'\Vw\ 2 + d'\ w afi ifrf. 

Suppose that the minimum of the second derivatives of w is taken along a direction a at (0, to). 
Then, is an eigen direction of the symmetric matrix D 2 iv(0, to). Thus, at (0, to), we have 

Z = IVaa + ef(t ) = 0, 

w^p = iiaj^p and ep e R", 
D 2 Z > 0, AZ > 0, VZ = and Z t < 0. 

In addition, we get 

C~[Ci ~h C ' ) 

(W 2 ) = 2ncl + — - and V x w ■ V x \i]\ 2 = 2c H w H at (x , f ). 
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Since au is uniformly positive, by (|2.7|l , we have 

-ef >2pa l i-w m jw m - ^pchx aa w m - eipa l i (|rj| 2 ) - 2x> s p<ihoijU> m 

- pafajeip (\i]\ 2 ) k + pdhvijW^ + p'a'hvjjW-^ + 2p , w d ^.w i fi> m 

+ 2/?|Vw| p " 2 (fe'zifr - C„b) W^as + 2 (c'wg - C«c) W^W^ 

- eip \pb\Vwf~ 2 + 2c] (Vw • V|t]| 2 ) + [pfo|Vwr 2 + 2c] a^L + p(p - 2)&|Vttf|P- 4 w 2 -wjL 
[b'lVwIP + c'IVwp + d'Jw^ 



(2.8) 



+ 



:= - Keip - Le 2 il> 2 



where 



K-i =2paJ.w mj - 4pcV a - 2Capa*Wij + p'a l 'wij + 2p'w^dl.Wij 

+ 2p\Vivf~ 2 (b'w-a - c-ab) w-a + 2 (c'lv-a - c^c) w w + b'\Vwf + c'\Vzv\ 2 + d' 

+ $ ^2ncl + C " (C ' 2 +C/) | + 20^ [pb\Vwr 2 + 2c] 



and 



L = -pdLwij - [pb\Vwf~ 2 + 2c) - p(p - 2)b|Vz«|P- 4 w|. 

By the boundary condition, the minimum point 0, to appears far away from the boundary. Hence, 
we can choose a compact subset S in Q such that 

(xo(f)/ zo(£)) e S and w(x, t) : bounded above and below in S. 

The equation (|2.1|) , when restricted on S, is nondegenerate. Therefore, the classical estimates for 
linear parabolic equations give 

\Dw(0, t)\ < C|M| L » (S ) and \D 2 iv(0, t)\ < C\\iv\\ L ~ {S) 

for constant C > 0. Hence, the quantity K and L are under control, i.e., there exists a constant M 
such that 

\K\,\L\<M at(0,f ). 
Hence, if we choose the function having the property 

xp t {t) > \KVp + \L\xp 2 , 

then, contradiction arises. Therefore, the time t — to satisfying (I2.3I) doesn't exist. Thus, 

inf inf wrr > -e sup ib(s), Vf > 0. 

xeQ,se[0,f] ^eR'Me^l ^ < s < f 
Letting e — > 0, we can get a desired conclusion. □ 

3. Geometric Property of Gelfand Problem 

In previous section, we discussed the positivity of the second derivatives of solutions for de- 
generated parabolic equation. It is a very useful tool for investigating the geometric properties of 
solutions to Gelfand's problem. We now address the long-time geometrical properties of solutions 
for the initial value problem with exponetial growth: 

(3.1) {e u ) t = am + Ae" in Q = Q x (0, oo) 
posed in a strictly convex bounded domain QcR" with 

(3.2) u = on Q, u > din Q 



GEOMETRIC PROPERTIES OF GELFAND'S PROBLEMS WITH PARABOLIC APPROACH 



and inital data 

(3.3) u(x,0) = Uq(x) < <p(x) 

where (p is the minimal solution of the Gelfand's problem \GP, 



In this section, we aim at providing the /-convexity of the minimal solution (p to the problem 
{ GP\) which also satisfies the boundary condition fll.4) , i.e., we want to show that 

e -2<pM : convex with respect to space variables. 



If we try to show the /-convexity of (p in ( GP,\ \, we can put v — e if and replace (p by -2 log v in 
the equation. Then 

VAV - |Vw| 2 - - = 0. 

2 

By direct computation, we get 

vAiv aa + 2v a Av a - 2Vv ■ Vv aa - 2Vv a Vv a = 0. 

Unfortunately, there are many terms, for example Av a , which are out of control. We don't even 
have any information about v aa . Hence, infering the geometric properties of (p from the equation 
( [GPD directly is very hard. However, according to the arguments in Lee and Vazquez's paper, 



| LVJ, the geometric properties of solutions to the nonlinear elliptic problem can be obtained by the 
geometric properties of the solution to the corresponding problem with parabolic flow. To apply 
their arguents to the Gelfand's problem, we start by showing the relation between the solution u 
of (33)-(33) and <p of §GPj\ . 

Lemma 3.1 (Approximation lemma). Let u(x, t) be a solution of (13.11) and let cp be the minimal solution 
of ( |GPa| ). Then, we have the following properties: For any sequence {t n }°° = -i with t n — » oo, we have a 
subsequence {t nk }T-, such that 

(3.4) lim \u{x, t nk ) - <p(x)\ -> 
uniformly in compact subset o/Q. 

Proof. Define the functional F(cp) by 

F(cp) = -f (IcpAcp + A^^dx 
and let g(t) = F(u(-, t)). Then a simple computation yields 

(3.5) g'(t) = -A [ e u{x ' l) [u t (x, t)] 2 dx < 0. 

Uniformly ellipticity of the cefficients in d3.ll > show that L e"^) dx is bounded for all t > 0. Hence, 
for some constant M > 0, 

(3.6) g(t) >-M Vt> 0. 

Therefore, by (13.51 1 and (13.61) , lim t ^oo g(t) exists and g'(t) — » 0. Hence, for any sequence of times 
(M~ =1 , tn -» oo, g'(t n )^0. 

Observe that the equation (13.11) is uniformly parabolic in Q. Thus, by the comparison principle, 
there exists a uniform constant C such that 

(3.7) < u(x, t) < <p(x) < C. 

Moreover, by the Schauder estimates for parabolic equation I1LSUL the sequence «(•, t n ) is equi- 
Holder continuous on every compact subset K of Q. Hence, by the Ascoli Theorem, the sequence 
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{«(-, t n )) has a subsequence {u(-, t nk )} that converges to some function h uniformly on every compact 
subset of K which is non-trivial. The readers can easily check the non-triviality of h. 
Multiplying equation (|3.1|) by any test function i] e C^°(Q) and integrate in space. Then for each 

tn 



(3.8) A I e u u t r\dx- I uAndx + A I e"r]dx. 

Jk Jk Jk 

Since the absolute value of the left hand side of 03.8I ) is bounded above by 



i 2 



(3.9) A (J* ife u dxJ^J e u (u f ) 2 dx^' 

and the second term of (13.9I ) has limit zero as t„ k — > oo, we get in the limit t, lk — > oo 

0= I hArjdx + A I e h r\dx, 
Jk Jk 

which is weak formulation of the equation 

(3.10) A/z + Ae'' = inK. 

By the arbitrary choice of a compact subset K in Q, ( 13.101 ) holds in Q. Therefore, h is a weak solution 
of the Gelfand's problem which is smaller than the minimal solution <p because of (I3.7D . Since (p is 
the minimal solution of Gelfand's problem, u(x, t„ k ) converges uniformly on every compact subset 
of Q to (p as t nk — » co and the lemma follows. □ 

3.1. geometric property. We first establish some estimates for the solution u of d3.ll ! which plays 
an important role for the geometric properties of solution u on the boundary. 

Lemma 3.2. Let u e C 2 (Q) x (0, oo)) be a solution of fl3.lH3.3l ). Then, there exists a constant C > such 
that 

\u TkV (x, t)\ < C\u v (x, 01, Vx e dCl (k = 1, • • • ,n - 1) 
where v and are normal and tangential directions to dD, at (x, t), respectively. 

Proof. Without loss of generality, we may assume that x = and outer normal direction e v = e n . 
For any 1 < k < n — 1, we consider the directional derivative 

d Tk u = P(x) • Vw, (P(x) = ( Pl (x), • • • ,p n (x))) 

which is the same as a tangential derivative on the boundary <9Q. By the strictly convexity of 
domain Q, we can find the largest circle of radius R — that touches the domain Q from inside 
at x = in the plane generated by two directions e %k and e v . Note that the reciprocal of R lies on 
between principle curvatures at x = 0. Moreover, (R - x n ) + x^u n is a tangential derivative on 
the circle. Thus, the directional derivative dj k u can be expressed in the form 

d Tk u{x, t) = (R- x n + 0(|x| 2 )) u k (x, t) + (x k + 0(|x| 2 )) u n {x, t) 

near the point x = 0. 

For a constant < p < 1 and a time to > 0, let B p = B p (0) be the ball of radius p centered at and 
define the new functions h+ and H+ by 

2 

h+ = dj.u + > (dj.u) , 

(3.11) " * i~t { 1 

H+ = r\ 2 h+ + px 2 , 
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where rj(x, t) - (p 2 - |x| 2 ) (t - to + p 2 ). Then, it can be easily checked that 

2 

(3.12) \ H± ( X/t )\<-EL Xn _ 

on {e?(Q n Bp)} x [t - p 2 , t ] and (q n B p ) x {f = t - p 2 } 

Let us now show that LH+ > and LH_ < in Q n Bp where L is defined by 

Lf = Af + Ae u f-f t . 

By direct computation, we can find a constant pi > such that 



n-l 



n-l 



£ \d Tl u\ 2 > y L l Vu ^l 2 - Ci (l + |Vu| 2 ) VO < p < Pl 
/=i /=i 
for some constant C\ > 0. Thus, we have, for some constant c > 0, 

LH + > ;] 2 (a/i + + Ae !i /i + - (h + ) t ) + 2h + (ryAr? + |Vr?| 2 - rjrjt) 
+ 477V77 • V7/+ + 2p 

>2 n 2 ^d Tl u\ 2 
1=1 

[ n-l ( ( n n-l 

- crAl + 2 (R|V^«| + p\u n \) Yj \ l V "l + P Tj |mkI + Zi + |Ax ' M| + Ae " + ,Mfl 

- Arfe u (R 2 \V x <u\ 2 + 2Rp|V^w||w„| + p 2 |u„|) 

- 4n(n + 3)p 6 (R|V x ,u| + p\u n \ + R 2 \V x <u\ 2 + 2Rp|V^w||w n | + p 2 |w n | 2 ) 

n-l 

- 16p 3 77(R|V. T , U | + p\u„\)Yj |v(^t,")| 



- 8p 4 ?] 



n-l n-l 

Y \u ni \ + |A^m| + Ae" + \u t \ - 8Rp 3 t] Y |uju| - 8p 3 r]|Vu| + 2p 

(=1 !=1 



> 



n-l n-l 
if \VdT,u\ 2 - 16p 3 r] (R|V^U| + p\u n \) Y |v (d Tl u)\ + | 
/=l /=1 



+ 



+ 



Yu l Vu 'l 2 - 3 W 1 + 2CR|V*tt| + pKI) £ |Vu,| + £ 
Z=l /=i 



2_2 «-l 



£ |Vw/| 2 - 8p 3 rj (R + p)Y I V "/l + f 
1=1 1=1 

+ I - 2p 6 [Ko + K 1P + K 2 p 2 + K 3 p 3 + R 4 p 4 ] 



where 



K = 4n(n + 3) (R|V*/w| + R 2 |V x ,w| 2 ) 
Kt = 4 [n(n + 3) (|u„| + 2R|V y ui|w n |) + 2|V«|] , 

X 2 = y (l + |Vu| 2 ) + 4n(n + 3)|u fi | 2 + AR 2 e u \V x >u\ + 8 (Ae 1 ' + \u t \) + c(n - 1)(1 + 2R|V*>«|)|Vj 
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K 3 = c(n - 1) [|tt w ||V«| + (1 + 2R\V x >u\) (Ae u + \u t \)] + 2RApe u \V x >u\\u n \ 

and 

X 4 = c(n - l)\u„\ (Ae u + \u t \) + Ae"\u n \. 
We now choose a constant p2 such that , V0 < p < p2, 



p|Vw| 2 < 1, 2 7 (R\V x ,u\ + p\u n \f < |u»(0, t )\, 



36c 2 p 4 



R 2 



,2 



and 



l + 2(R|V^ M | + pK|)] <\u n {Q,t Q )\, 28 ( 1 + f) ^ 



K + K lP + K 2 p 2 + K 3 p 3 + K 4 p 4 < 



in jo Pi Bp} X [fo - p 2 , fo]- Thus, for < p < min {pi, pz], 

(3.13) L(H + )>0 ifp>max{K(0,f )lp 6 ,2 9 p 6 } = 2Ci|w„(0,t )|p 6 . 

Hence, by d3l2l and (l3l3l . 

H+(x, < -CxKCa ^o)lp 8 ^ on Q n Bp x (f„ - p 2 , £ ]- 
This immediately implies that 

h+(x, t) < C\U n {Q, t )x n in Q n Bp x |t - p fo 

Similarly, we can also show that 

h-(x, t) > -C 2 u n (0, t )x„ inQnBp x|t -p t 

for some constant C2 > 0. Therefore 

-CRu n (0, to)x n < V-(x, t) < v+(x, t) < CRu n {0, to)x n in Q n Bp X ^0 - X' f ° 

for constant C = ^ max{Ci, C2}. Taking the normal derivative e? n at (0, fo), we obtain 

(f+)»(<Uo) 



l«ibf(0/ fo)l 
and the lemma follows. 



R 



< C\u n (0)\, VI < k < n - 1 



□ 



Next, coming to our subject, we have the following result about preservation of f- convexity, 
which is easy but allows to present the basic technique. Our geometrical results will be derived 
under the extra assumption that Q is strictly convex. 



Lemma 3.3. Let Clbea strictly convex bounded subset in R" and assume that u is a solution of (|3.H> - (|3.3 
with the boundary condition \1A) . Let u = -2 log w. Then, for every t > 0, as x — > xq e dQ 

(3.14) w aa {x, t) = ^e"3" {^u 2 x - u aa ) >5 >0, (a = e a e R M , and \e a \ = 1) 

for a uniform constant 5q depending on <9Q. 
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Proof. By direct computation, we have 



n 



IV r. 



i(x, t) = J" (-ul - u aa ) , on dQ. 



If e a = x, a tangent direction at xq to dQ, then u a = 0. Hence, we need to estimate u aa . For this, we 
use the fact that dQ is strictly convex. Without loss of generality, we may assume that Xq - and 
the tangent plane is x n — 0. We may also assume that the boundary is given locally by the equation 
x n = f(x'), and x' = (x\, • • • , x n -\). We introduce the change of variables 

(3.15) x a = x a , x n = x n - f{x'), g{x', Xn, t) = u(x', x n , t). 
Then along tangent directions we have 

W-aaipC / %n, = £aa(x , %n, — ^-§ax„ , %n, fyfa(x ) 

+ gx„x„ (x',X n , t){f a {x')) 2 - gg n (x', X n , t)f aa {x'). 

Since f a (0) = and f aa (0) are nonzeros along all tangent directions, we get 

(3.16) u aa (0, 0, t) = -g Xn (0, 0, t) f aa (0) = -u Xn (0, 0, f) f aa (0). 

Since u is a solution to a parabolic equation with uniformly elliptic coefficients, by the Hopf's 
lemma for the classical partial differential euqation, < cq < |Vu|. Hence 

u Tt (x Q/ t) = -u v (x 0/ t)f TZ (0) < 

and 

Co/tt(0) 



(3.17) w TT (x, t) = hi 5" [-u 2 . - W TT J 

Let e a - v. On dQ 



> 



as X — > Xq 6 <9Q. 



= A (e") f = Aw + Ae" = am + A. 



Thus, we have 



- w vv = + A + / i u TiTi = -M^ + A + (n - 1)m v H(<9Q) on <9Q, 



!=1 

where H(dQ) is the mean curvature of <9Q at Xo- By the boundary condition (J1.4|) , we can obtain 
(3.18) w vv {x, t) = -e^ u [-ul - u w ) = - (-u* + A + (n - l)u v H(<9Q)) > on dQ. 

Finally, we check the case that the minimum of the second derivatives of w occurs along a 
general direction e a at xq. Since the outer normal direction e v is vertical to the tangent plane, we 
can express a general direction e a by 



(3.19) 



e a = he aT + k 2 e v , [k 2 + k\ = l) . 



where e« T is a direction contained on the tangent plane to the graph of u at (0, t). Hence, the second 
derivatives of w can be written in the form 



iv n 



2 

Thus, at x — xq, we have 
(3.20) w a 



k 2 
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By Lemma [3^21 there exists some costant C\ > such that 

(3.21) \u arV (0, t)\ = |zvv(0, f)| < Q \v v (0, t)\ = C 2 \u v (0, t)\ 

Combininig ( 13.211) with d3.20|) , we can get 



k\ f-w 2 , + A + (n - \)H{dQ)u v \ - dhk 2 \u v \ + k\f ax a T u v 



2 

Therefore, if u satisfies the boundary condition \1A\ for 



at x = Xq. 



4|/« T « T (0)|' 

then the result (13.141) holds for all e a e R", (\e a \ - 1) and the lemma follows. □ 

Theorem 3.4. Let Q be a convex bounded domain and let Uq > be a continuous and bounded initial 
function which satisfies (|3.3|) . Then, the solution u of (|3.1|) -(|3.3 | ) with the boundary condition (|1.4|) z's 
f -convex in the space variable for all t > 0, i.e., 

(3.22) D 2 (e-5 U )>0. 

Proof. Let if = e~? u Then, the new function w satisfies 

(3.23) Awt - w Aw - w\Vw | - — w. 

Hence, w is a solution to the equation (|2.1|) with p(w), fl, ; (Vw), b(w), c(w) and c?(tt> ) being replaced 
by w 2 , I n , 0, —if and -jW respectively. Here, l„ is the nxn identity matrix. 

We now are going to check that the condition 1.1 and 1.2 given in Section|2]hold for the coefficients 
in the equation ( I3.23D . It is trivial that 

b(w) = 0, c(w) - -w, d(w) - -^w : convex. 
In addition, by direct computation, it can be easily shown that 

P 2p 2w 2 

Hence, the equation (13.231 ) have the condition I and II. On the other hand, Lemma 13.31 tells us 
that the solution of the equation (13.231 ) with zero boundary condition is convex on the boundary. 
Therefore, by Lemma l2Tl it is also convex in the interior of domain Q and the lemma follows. □ 

Corollary 3.5. IfQ. is convex, the the stationary profile <p(x) of u{x, t) is f -convex, i.e., 



Proof. Take the inital data as before. By the asymptotic result, Lemma 13.11 we have uniform 
convergence between <p(x) and u(x, t). Hence, the conclusion follows. □ 

Note that <p satisfies the equation 

A(p + Ae v = in Q. 

Then w = f(<p) 2 = e~f satisfies 

— — i —i 2 — 

iv Aw - |Vw| - Azv = in Q. 

It has the similar form to the equation in (3.12) of the paper ILVj . Hence, following the same 
arguments as in the proof of Lemma 3.6 in IILVL we get an improved result. 
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Lemma 3.6 (Strictly /-convexity). IfD. is smooth and strictly convex, then the minimal solution <p(x) is 
strictly f -convex: there exists a constant c\ > such that 

D 2 f((p) = D z (e-y)<-c 1 I. 
The constant c\ depends only on the shape ofQ. 

4. Boundary Condition for the Geometric Properties 

Through the previous sections, we investigated the /-convexity of the minimal solution of 
Gelfand's problem, ([GPa), under the assumption that the minimal solution has the property (|1.4|) 
which we call boundary condition. As mentioned before, if we have the boundary condition 
removed from our setting, it is very difficult to get the /-convexity of the minimal solution of jGP^ 



on the boundary. Thus, we had no choice but to add the boundary condition (|1.4[) for the result. 
However, the problem is that we couldn't guarantee the existence of the minimal solution of ( |GPa| ) 
having the boundary condition {LA) . Therefore, we need to check whether it is possible to satisfy 
conditions, ( |GPa| > and (|1.4|) , simultaneously. Otherwise, it is meaningless to apply our geometric 
results in physical models. 

Before we finish this work, we will introduce some properties related the boundary condition 
(|1.4|) in this last section of the paper. We are now ready to state and prove the main result in this 
section. 

Lemma 4.1. Let CI be a ball with radius r > 0, i.e., Q = B r , and let ip,\ and be solutions of the Gelfand's 
problem \GP\) in a ball B r . 

(1) Suppose that > (p,\. Then, there exists a constant ro > such that ifG{§\, A, B r ) > then 

G(xp A ,A,B r )>0 Vr>r . 

(2) Suppose that ip,\ > cp A . Then, there exists a constant r\ > such that ifGiip^, A, £>,•) < then 

G((j) A ,A,B r )<0 Vr>n. 

(3) For any < A < A*, there exist a constant r2(A) = r^ > such that 
(4.1) G(^ A/ A / B r )>0 Vr>r 2 . 

(4) For any < A < A*, there exist a constant r^A) - r$ > and a solution x^x of the Gelfand's problem 
defined on B r3 such that G{\^\, A, B r3 ) > fails. 

Proof. For the convenient, we assume in this proof that B r = B r (0) a ball of radius r which is centered 
at 0. Define the function 9,\ : R" — » Rby Q\ (x) = ^ (r 2 - |x| 2 ). Then, by a direct compuation 

A6(x) = -A. 

For any solution (p\ of the Gelfand's problem ( GP\) , by the maximum principle for the super- 



harmonic function, we have (p\ > 0a in B r and <p\ - 6\ on dB r . This immediately implies that 

Ar 

(4.2) (px,v<— ondB r . 

(1) By assumption, we have 

^)\>(\)\> 6\ in B r and ty\ - (p,\ - 9\ on dB r . 
Thus, it follows that 

Ar 

ty\,v < (pA,v < on dB r . 

n 
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where v is the outer normal vector at the boundary dB r . Moreover, it also follows from 
G((pA, A, B r ) > that the mean curvature H(B r ) satisfies 

1 \ 

(n - l)H(dB r ) + K< --4> A/V - — on dB r . 

2 <pA,v 

Thus, for r>^r-/ 

Gty K , A, B r ) =\ ($x,vf + A + (n - l)i/> A ,vH(<9Q) + Ki[> Aa , > i^ A , v (i/> A , v - <£ A/V ) + A |l - 

Hence, G (i^ A/ A, B r ) > holds for r = ^if. 
(2) In this case, the mean curvature H(B r ) satisfies 

1 A 

(n - l)H0B r ) + K > --ip A/V - — - on dB r . 

2 </> A/ v 



Thus, G ^ A/ A, B,) < also holds for r\ = J ^j- since 
G(c/> A , A, B r ) =i (<p A , v f + A + (n - l)cj> A , v H(dQ) + K<£ A/V < hp A , v (</> A , v - i/> A , v ) + A (l - ^ j 
^ - ^(i^ ¥ - A) < A (*M^(£ - lU o. 



i/>A,v /\2 ' ' / \ A/V /\2n 2 
(3) Suppose that \p A is a solution of the Gelfand's problem. Then, it follows from (|4.2[) that 

Ar 

u»A,v ^ on o>B,.. 

n 



Therefore 

1 1 T Ar 91m -11 

+ A 



G(i/> A/ A, B,.) = ^ A/V [t^ A ,v + 2(n - l)H(B r ) + K] + A > ^a,v 



n r 



> iir >^V + f. 

Hence, 63) holds for r 2 = + f . 

(4) By a simple computation, it is easily checked that the function g(x) - log y ^a ) is a solution 
of Gelfand's problem on B i . In addition, we also have 

V-i- 

1 {§v f + A + (n - l)g v H(B r ) + Kg v = + A - + K Sv < on <9B 

Hence, 



This finishes the proof of (4) and the lemma follows. 

□ 
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We give next the generalization of the previous lemma for a strictly convex bounded domain Q 
with smooth boundary. For x e <9Q, we denote by 

Tq(x) - sup{r e R : B r c Q and £>,. n dD, - {x}} and tq = inf tq(x). 

xedd 



Then, we easily conclude that 



and 



H(dO) < — 



\p\ < on dll 



for any solution \p,\ of the Gelfand's problem |GPa). Hence, following the similar computation as 
in the proof of Lemma |4~T1 we can easily extend the results for the boundary condition on a ball to 
more general boundary setting. 

We finish this work with stating the following result. 

Theorem 4.2. Let CI be a strictly convex bounded subset of R" and let \p\ and <p A be solutions of the 
Gelfand's problem ( |GPa| > in Q. 

(1) Suppose that ip,\ > tp^. Then, there exists a constant r$ > such that ifG(<p Ar A, Q) > then 

G(i/>a,A,Q)>0 Vr Q >r . 

(2) Suppose that ip,\ > tp^. Then, there exists a constant r\ > such that ifG(\p\, A, Q) < then 

G(^ A/ A,Q)<0 Vr Q >r 1 . 

(3) For any < A < A*, there exist a constant r2> such that 

G(xp A ,A,Cl)>0 Vr n >r 2 . 
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